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A one-dimensional  theory  based  upon  fundamental  flow  relationships 
is  presented  for  analyzing  the  behavior  of  one  or  more  gas  streams  flowing 
through  a single  nozzle.  This  conpound-corpressible  flow  theory  shows  that 
the  behavior  of  each  stream  is  influenced  by  the  presence  of  the  other 
streams.  The  theory  also  shows  that  the  behavior  of  conpound-compressible 
flow  is  predicted  by  determining  how  changing  conditions  at  the  nozzle  exit 
plane  affect  conditions  within  the  nozzle.  It  is  found  that,  when  choking 
of  the  corpound-conpressible  flow  nozzle  occurs,  an  interesting  phenomenon 
exists.  The  conpound-corpressible  flow  is  shown  to  be  choked  at  the  nozzle 
throat,  although  the  individual  stream  Mach  numbers  there  are  not  equal  to 
one.  This  phenomenon  is  verified  by  a wave  analysis  which  shows  that,  when 
choking  occurs,  a pressure  wave  cannot  be  propagated  upstream  to  the  nozzle 
throat  even  though  some  of  the  individual  streams  have  Mach  numbers  less 
than  one.  Algebraic  methods  based  on  this  compound-compressible  flow  theory 
are  used  to  demonstrate  the  usefulness  of  this  approach  in  conputing  the 
behavior  of  compound-compressible  flow  nozzles.  A comparison  of  the 
compound-compr^  sible  flow  theory  with  three-dimensional  computer  calculations 
shows  that  the  effects  of  streamline  curvature  on  nozzle  behavior  can  be 
disregarded  for  many  practical  nozzle  configurations.  Test  results  from  a 
typical  two-flow  nozzle  show  excellent  agreement  with  the  predictions  from 
the  theory. 
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COMPOUND-COMPRESS IBLE  NOZZLE  FLOW 


Introduction 

Modeun  propulsion  engines  often  exhaust  several 
different  streams  of  gas  side-by-side  through  a single  nozzle,  Fig. 
1.  These  flows  ran  exhibit  sizeable  eompressiblity  effects  and 
they  will  la*  referred  to  here  as  eoni|auuid-eompressible  nozzle 
flows.  The  pur|x»se  of  this  pajier  is  to  provide,  for  the  first  time, 
a simple  method  to  prediet  the  l>eh!ivior  and  clarify  the  under- 
standing of  such  flows. 

A uue-dimensioual  analysis  similar  to  that  used  in  single-stream 
compressible  flow  problems  is  applied  here  to  com|>ou nd-com- 
pressible  flow  problems.  The  great  advantage  of  this  type  of  ap- 
proach i>  that  it  provides  physical  insight  into  the  nature  of  the 
flow.  1 

Mixing  lietwmi  the  various  streams  is  not  considered  in  the 
development  of  the  bade  theory,  but  its  effect  on  eontjsnind- 
compressible  flow  U*liavi<»r  will  la*  discussed.  It  will  Ik?  shown 
that  mixing  often  hits  u negligil  le effect  on  the  flow  lieliavior. 

The  usefulness  of  the  ront|>uimd-coiiipressihie  flow  theory  is 
demonstrated  by  comparing  its  prediction*  with  Imth  three- 
dimensional  computer  calculation*  and  ex|K*riiuciital  results. 

'Hie  basic  approach  u«**d  to  develop  the  coiu|s>imd-eompreasi- 
ble  flow  theory  will  lie  to  determine  how  changing  conditions  at 
the  nozzle  exit  plane  change  conditions  within  the  nozzle.  This 
will  Ik*  seen  to  lie  tlie  heart  of  the  matter  and  all  results  obtained 
m this  pa|K*r  are  presented  in  tin*  light.  Note  should  la*  taken 
of  a pioneering  contribution  to  «'oni|*iuml-4s  impressible  nozzle 
flow*  made  by  IVarson,  Holliday,  ami  Smith  |l}.*  Their  results 
are  consistent  with  the  general  eoneliisioits  arrived  at  in  this 
paper. 
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OiiOiamional  Coap«nJ-Coaprtuikli 
Nuzli  Flew  Theory 

The  development  of  one-dimensional  compound-compressible 
nozzle  flow  theory  follows  that  of  Shapiro  [21  for  single-stream 
flow.  The  most  important  alteration  is  that  the  fluid  statu* 
pressure  is  chosen  as  the  dc|tendent  parameter  lierause  it  can  vary 
only  along  the  nozzle  in  one-dimensional  flow,  w lie  re  as  all  other 
fluid  properties  can  also  change  from  strram-to-stream  acrots  the 
nozzle. 

This  analyst*  is  sufficiently  general  to  include  any  arbitrary 
mimlier  of  streams  designated  by  the  integer  n.  For  examplp, 
at  any  position  in  the  nozzle, 
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where  .1  is  tlie  total  flow  aira.  .1 , is  the  flow  area  of  tlir  it h stieain, 
and  x is  tiw  axial  nozzle  |s>»ilmn  csiofdmate.  In  single-'tieam 
<:ne~d;mcu-tniini  llieory,  »M  if/  is  arbitrarily  small  ami  this  earries 
oxer  into  the  present  rase  whrte  all  il.t,  if/  are  arbitrarily  »nudl 
The  l ran* verse  pressure  gradients  caused  by  streamline  nirvati.rr 
ran  then  I**  neglected  and  thi*  lead*  to  tlie  ronrbtsmn  that  static 
pressure  uoniy  a function  of  axial  |*:sitmn. 

Jt  i»  also  assumed  that  the  flow  in  earh  stream  is  *tra»|v,  adia- 
itatie,  and  tsrntmpie  and  that  each  fluid  i»  a perfect  gas  with  con- 
stant thermodynamic  pmprrtirs.  Note  that  thesr  assumption* 
exclude  mixing  effect  s. 
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to  Fife.  -,  the  consequences  of  the  choice  of  inlet  pressure  can  now 
he  examined. 

If  the  inlet  pressure  is  sufficiently  high  (curves  a and  b),  equa- 
tion (6)  will  yield  values  of  M<*  small  enough  that  0 will  be  >0 
everywhere  in  the  nozzle  and  p will  therefore  change  in  the  same 
direction  as  A throughout,  equation  (3).  In  particular,  both  p 
and  A will  have  their  smallest  values  at  a geometric  throat  where 
A reaches  its  minimum.  Note  that  the  integration  also  shows 
what  the  back  pressures  must  be  to  maintain  these  flows. 

At  the  same  time,  the  differentiation  of  0 with  respect  to  p 
yields 


d0 

dp 


Aj 


(1  - M,»)»  + 2 


kj  + 1 
2 


> 0 
(7) 


Pig.  2 Cempmmd-cemprestlM*  tow  in  • nasxt*  at  tow!  f«w«try 


Each  stream  may  then  be  separately  treated  as  a single-stream 
one-dimensional  flow  (Shapiro  (21,  Table  H.2).  Consequently, 


where  it,  is  the  ratio  of  s|>enfir  heats,  M,  is  the  Mach  number  (M, 
- VJy/klRJ'i)  for  the  t'th  stream,  and  p is  the  fluid  static  pres- 
sure. Equations  ( ! ) and  (2)  may  be  combined  to  yield 


where 


(3) 


(4) 


The  term  0,  the  compound-flow  indicator,  will  *iil*eqiieuily  lie 
shown  to  be  significant  in  determining  the  nature  of  the  flow  (i.e., 
whether  it  is  rompound-subaonir  or  row pound-au peraon ir). 

The  behavior  of  compound-comprrssiblc  flow  in  a uoxxUs  of  fixed 
geometry  (an  lie  most  profitably  examined  by  regarding  the  inlet 
pressure  as  an  independent  variable,  Fig.  2.  If  the  stagnation 
pressures  (pa,),  the  stagnation  temperatures  (7*#,),  and  the  gas 
properties  are  constant  and  known  for  each  stream,  the  mass 
flows  (tra)  are  functions  only  of  (lie  local  pressure  and  the  local 
flow  area: 


A.p*. 

>/T 


‘Tins,  for  any  given  value  of  p at  the  inlet  plane,  where  the  .4* 
are  known,  equation  <’»  j may  l*e  used  to  determine  the  corre- 
sponding values  of  r,.  With  the  ir,  fixed  by  the  inlet  prewar*, 
it  ran  la*  *een  from  equation  ».’»>  that  the  local  At  are  functions 
only  of  the  local  p and  known  quantities  in  the  remainder  of  the 
nubile.  Since  the  local  M.  are  also  functions  only  of  the  local  p 
and  known  quantities,  namely. 


it  follows  direr*  !>  that  I hr  local  value  of  £#,  equation  (4).  is  a 
function  only  of  the  inlet  pressure,  the  local  prewiire,  and  known 
quantities  Therefore,  equals  si  i.l)  ran  f»  integrated  in  principle 
from  inlet  to  exit  for  any  rhosm  value  of  inlet  pressure.  ttrfemng 


which  shows  that  0 always  changes  in  the  same  direction  as  p. 
Therefore,  for  curves  a and  b,  (3  will  also  change  in  the  same 
direction  as  .4  and  will  also  have  its  minimum  value  at  the  throat. 

As  the  inlet  pressure  is  decreased,  the  value  of  (3  at  the  throat 
will  also  decrease.  In  fact,  when  the  inlet  pressure  is  chosen  to 
be  sufficiently  small,  0 reaches  zero  at  the  throat.  When  this 
occurs,  equation  (3)  is  indeterminant  and  no  longer  serves  to 
determine  the  axial  pressure  gradient  at  the  throat.  Under  this 
condition,  application  of  L’ Hospital's  rule  to  equation  (3)  yields 

l (ln  p) 


(8) 


The  geometry  of  any  throat  is  such  that  rfM/dr*  i*  always  >0. 
Therefore,  rf{ln  p)/<ir  will  be  either  the  positive  or  negative  root 
of  a real  number. 

Curve  c represents  the  choice  of  the  positive  root  while  curve  d 
represents  the  choice  of  the  negative  root.  Comparison  of 
curves  c and  d reveals  a familiar  single-stream  compressible  flow 
situation:  The  geometric  throat  is  a saddle  point  for  two  isen- 
tropir  solutions  in  the  divergent  section  of  the  nozzle.  No  bark 
pressure  between  curves  c and  d ran  correspond  to  an  isentropie 
flow.  It  is  anticipated  that  those  hack  pressures  which  do  not 
correspond  to  iarntropir  solutions,  such  as  that  of  curves  may 
lie  reached  by  means  of  compound  akork*  initiated  at  some  point 
on  curved. 

The  hehavior  of  the  flow  along  curve  c is  similar  to  that  of 
curves  a and  h.  i.e..  A,  p,  and  0 will  reach  their  minimum  values 
at  the  throat.  Note  that  the  positive  root  of  equation  (8)  is 
chosen  only  when  the  liaek  premure  corresponds  exactly  to  that 
of  curve  e. 

The  implications  of  the  choice  of  the  negative  mm  of  equation 
(K)  will  now  he  considered  in  detail  (curve  d ).  Hincr  equation  (7 ) 
has  shown  that  $ always  changes  in  the  same  direction  a»  p,  d0/dt 
must  also  be  negative  at  the  throat.  Accordingly.  0 must  de- 
crease from  positive  to  negative  as  it  passe*  through  seen  at  the 
throat.  Furthermore,  with  0 negative  entering  the  divtrgwit 
section  i where  d.t  ids  is  > tij.  simultaneous  examination  of  equa- 
tions <:t>  and  1 7)  »hows  that  the  local  value*  of  both  p and  0 must 
continue  to  tlrrreasr  through  the  divergent  section.  Note  that 
the  integration  of  eqiutfkor.  still  shows  what  the  bark  pressure 
must  1st  l»»  maintain  thi»  flow. 

No  iwmtroptr  »JuIss»k  exist  for  snlH  pnwnuw  rurrespsnding 
U»  value*  le**  than  that  of  curve  d lectiue  0 would  reach  arm 
upstream  of  the  geometric  thmat.  This  would  wswlt  in  an  infinite 
axial  pressure  gradient,  equation  (3). 
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Stum*  interesting  conclusions  etui  Im?  made  by  examining  the  in- 
H uen re  of  bark  pressure,  pm,  on  the  inlet  pressure.  A backpres- 
sure grenler  than  that  of  curve  e will  atTert  the  pressure  at  the 
inlet  plane  and  will  thus  influence  the  flow  rates  of  the  individual 
streams.  Any  bark  pressure  less  than  that  or  curve  e will  affect 
neither  the  inlet  pressure  nor  the  flow  rates.  This  condition  wil) 
l>e  referred  to  as  compound-choking.  Under  such  conditions,  the 
nozzle  geometric  throat  controls  the  t>ehavior  of  the  flow. 

Also,  since  dp /dx  is  always  <0  for  curve  d,  Bernoulli’s  equa- 
tion shows  that  a ron'inuous  acceleration  of  the  flow  takes  place 
throughout  the  nozzle.  Applying  this  to  equation  (3),  it  can  be 
seen  that,  for  continuous  acceleration  of  the  flow,  > 0 whenever 
dA/dr  < 0 and  (i  < 0 whenever  dA/dx  > 0.  Thus,  for  con- 
tinuous acceleration  of  the  flow  in  a single-nXre&m  convergent- 
divergent  nozzle,  examination  of  0 reveals  that  the  flow  must  be 
subsonic  in  the  convergent  section,  sonic  at  the  throat,  and 
supersonic  in  the  divergent  section,  equation  (4).  In  the  follow- 
ing section,  it  will  lie  shown  that  for  compound-compressible  flow 
an  analogous  situation  exists:  The  flow  must  be  compound-sub- 
sonic in  the  convergent  section,  compound  sonic  at  the  throat, 
and  comjsnmd-supersonic  in  the  divergent  section.  It  will  also 
Ik*  shown  that  these  regimes  are  differentiated  by  the  com- 
pound-fluw  indicator,  0. 


Compound  Waves 

The  cuinponiul-chokiiig  phenomena  just  desrrilx'd  can  Ik*  ex- 
plained by  examining  the  effects  of  small  pressure  disturbances 
on  the  flow.  A diagram  of  micIi  u disturbance  is  shown  in  Fig.  3. 
ll  is  consistent  with  one-dimensional  theory  to  take  the  flow  area 
as  constant  in  wave  calculations.  If  a weak  plane  pressure  dis- 
turbance is  imimsed  on  the  flow,  Fig.  3(*j),  this  disturbance  cannot 
propagate  at  different  absolute  velocities  in  each  stream  without 
violating  the  condition  that  the  static  pressures  at  the  stream 
interfaces  lie  equal.  Therefore,  the  wave  must  Ik*  continuous 
and  must  travel  as  a single  ouu|kiuimI  wave,  Fig.  3(fc).  Although 
the  wave  is  not  ru*<*essarily  plane,  the  pressure  rise  neross  it  cannot 
vary  from  stream  to  stream. 

As  indicated  in  Fig.  3(fc),  the  absolute  terminal  velocity  in  the 
upstream  direction  of  the  compound  wave  is  designated  by  a.  It 
follows  directly  that'  a > 0 rormqmnds  to  eompound-autaonir 
flow;  a « 0 corresjHUids  to  .•ompound-aonic  flow;  a < (t  corre- 
s|M>nds  to  rompound-auperaonir  flow  . 

An  analytical  expression  for  the  compound  wave  velocity,  a, 
can  Ik*  derived  by  treating  each  stream  separately  as  a flexible 
lids-  and  conserving  mass,  momentum,  and  entropy  arross  the 
eom|SMiud  wave  in  the  frame  of  refereio'e  of  the  rom|snu>d  wave. 
It  follows  that 
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where  A signifies  the  change  acton*  tlie  com|«nintd  wave.  Since 
l 1m*  flow  area  is  constant, 
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3b  UPSTREAM  MOVING  COMPOUND  WAVE  IN  TERMINAL  STATE 
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Examination  of  equation  (12)  shows  that  the  compound  wuve 
velocity  a,  and  the  compound  flow  indicator  £#,  must  always  have 
the  same  sign.  Therefore  (j  > 0 corresponds  to  compntiud-sn)>- 
soiiie  flow;  0 <■  0 corresponds  to  eoni|M«nnd-sonie  flow;  $ < ll 
corresponds  to  compound-supersonic  flow. 

Xotr  that  comftound-choktng  can  onli y occur  at  the  nozzle  mini- 
mum area  for  it  is  there  only  that  (i  ran  equal  zero.  Note  also  that 
small-amplitude  compound  waves  cannot  move  upstream  in  th* 
compound-supersonic  region.  One  would  therefore  e\|K*et  that 
eoiii|MMind-shoek  waves  could  arise  and  cause  the  steady  rum* 
pound-compressible  flow  to  lie  uonisen tropic  in  the  c"Ui|>ntiiid- 
supersonic  region.  By  the  same  reasoning,  cuinpoii  ml  -shock 
waves  are  not  expected  in  the  eoin|Kiuiid-sul>soiiie  region. 

Several  conclusions  mav  Ik*  drawn  from  this  analysis: 

1 The  concept  of  flow*  choking  at  Mach  one  is  no  longer  \alid 
ill  compound-flow  analysis.  Indeed,  when  eoiit|xiutid  flow  is 
choked,  tlie  individual  stream  Mach  tnunliers  at  ihe  throat  will 
not  Ik*  equal  to  one  (except  for  the  unique  case  where  the  stagna- 
tion pressures  of  all  the  streams  an*  equal  and  tin*  Mach  ntimliers 
of  all  the  streams  are  unity  1.  Rather,  compound-choking  is  de- 
termined hy  the  compound-flow  indie itor  and  can  occur  only 
wlieu  d ~ 0 at  the  minimum  nozzle  flow  area. 

2 Not  ever*  tream  need  have  a Mach  numl>er  <1  in  order 
that  the  flow  Ik*  lompoiind-siibsoiiic  and  not  every  stream  need 
have  a Mach  numtier  '» i in  order  that  tin*  flow'  Ik*  rom|>omid- 
*uprr»ohir.  In  fart,  equation  (4)  show*  that  the  various  stream* 
influence  (i  in  pro|Mirtioii  to  their  flow  areas,  which  agrees  with 
intuitive  reasoning. 

3 Tlie  ronitHMind-flow  regimes  determined  by  d,  eunqiound- 
siilo.oi.ir  and  eoinja»un«l-so|ipr*.n»ic,  are  analogous  to  tlie  sul*- 
soiur  and  siqiersonir  flow  regimes  encountered  in  single-strram 
nozzles.  In  fact,  the  usual  tingle  -stream  results  arr  ofttained 
when  n — I in  the  foregoing  equations 


jr  s t.  - ia  - o ( io) 

across  l he  w a\v  And  »une  /»,  = p ami  Ap.  - fhp.  eqiiatarti*  (U) 
and  * l*li  max  Ik*  combined  to  vield  the  ilevred  relation  for 
tlie  rotlt|»»;nd  Wave  wlorilV  ft.  naWrlv, 


Equations  i 4 > and  ill)  mat  Ik*  is.mhinrd  to  yield 


Caapatatiaial  Pnatim  tor  Ciapul-Ciaptuikto 
Nintol 

1 W-sptlr  the  itml  f*»r  differential  equation*  in  drawing  rvmrlu* 
sion*  aUxit  ooe-dimrnsiotia!  compound-compressible  flows,  the 
*olut k>iis  for  ll»e  Iwhavior  of  these  flow*  will  fie  seen  to  require 
only  algebraic  computations 

Although  the  equation*  in  the  preceding  nrrtiotts  were  deuvr  1 
for  n flows,  to  illustrate  the  abdication  of  this  anah  »i».  only 
lWo-*trram  coiivergrtit-divrrgrul  nozzles  will  lw  rotisiderrd  m 
thi»  section.  *1111%  will  not  only  Itsmll  in  the  development  >4  less 
involved  equation*  but  Will  aha.  rtanfy  the  role  of  the  gesitnetrie 
throat  and  the  phenomenutt  uf  choking. 
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compound-  compressible 
flow  in  determined  by  the  relationships  of  k,,  I{„  .4  „ wtt  7'o„  pu„  and 
p.  Equations  ( 1 ) and  (5)  can  be  combined  to  yield 

- £ At  - A (13) 

i-t 

Using  equation  (13),  the  following  expression  may  be  written  for 
two  streams  at  any  point  in  the  nozzle: 
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These  are  the  fundamental  equations  for  solving  two-stream 

compound-compressible  flow  proliems. 

When  unrhoked  flow  exists,  it  has  been  shown  that  the  flow  ,se- 
havior  is  determined  by  the  back  pressure,  Fig.  2,  and  pm  - p.,,». 
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Thu*,  for  un choked  flow,  equation  ( 14)  may  be  applied  at  the  exit 
plane.  Typical  solutions  of  equation  (14)  for  unchoked  flow  are 
presented  in  uoudimensiotial  terms  in  Fig.  4 as  the  curved  lines 
of  different  /li*//1»it.  Although  these  curves  were  generated  for 
particular  values  of  k,  H,  and  temperature-corrected  mass  flow 
ratio  (v>t  y/T*/v>i  y/tQ i),  similar  solutions  can  lie  generated  for 
any  other  values  of  these  parameters. 

Under  choked  flow  conditions,  it  has  been  shown  that  the  flow 
liehavior  is  determined  by  the  nozzle  geometric  throat  where 
ft  - 0.  Equations  (4 ) and  (C)  can  lie  combined  to  yield  the  follow- 
ing equation  for  choked  conditions: 

u>»  Vtn  m (kt  Pw\  | *i  -J  T fp\  *'  _ .1  _ l 

«’i  VTo,  \*i  pw/  ( 2 LW  J ~ ) 

■i  -viar  •!; 


(16) 


Equations  (14)  and  (16)  may  lie  simultaneously  solved,  by  trial- 
Hiul-crr»r  determine  the  relationship  lietweeu  pn/p oi,  u> , y/f «/ 
ic,  y/U  * / .1  and  piwt/pm  at  choking  for  any  given  gas 

property.  Tl»c  relationship  between  pat/p D),  y/T*/ it,  y/T*, 

and  .4|*/.l»i.f.*t  is  shown  in  Fig.  It  can  lie  seen  that,  for  any 
given  value  of  u’:  y/T*/v\  y/f»,  there  is  a unique  pm/pw  corre- 
sponding to  each  .1 1 */ .1  ihriwi*  These  values  appear  as  horizontal 
lines  in  Fig.  4 liecniise  tlie  back  pressure  does  not  affect  tin* 
choked  solution. 

Referring  to  Fig.  4 for  any  given  combination  of  .1  »*/.!#*,»  and 
.1 1 */.4  in  a single  nozzle,  tlie  interaction  «»f  the  (<orres|a Hiding 
lines  indicates  the  onset  of  clinking  (curve  c in  Fig.  2).  For 
values  of  pm/p».  smaller  than  that  at  the  intersection,  the  flow 
liehavior  is  given  by  tlie  unchoked  curve.  For  values  of  p,,/p. 
greater  than  that  at  the  intersection,  the  flow  liehavior  is  given 
by  the  clinked  straight  line.  Thrrrforr,  the  entire  flow  behavior  of  a 
rompountl<omprrAMiUr  no;:lr  ran  be  dearnbrtl  by  a ainfflr  line. 
For  example,  the  dashed  line  in  Fig.  4 represents  a nozzle  with 
-t|*/.tik*.*i  » 0.431  and  .1  *'•/.!**.%  » 0.226.  I These  are  the  ac- 
tual dimensions  of  the  test  model  descrilied  in  the  exjierimentAl 
section  of  this  paper.) 

The  Uliavinr  of  M,  and  M-  at  tlie  throat  for  rhokeil  flow  ns  a 
function  of  it*  v/7‘w/wi  y/Tm  and  fixed  .ItV/tn,*.,  and  gas 
properties  is  shown  in  Fig.  6.  Sole  that  neither  at  ream  is  aonk  of 
the  nettle  throat. 
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Additionally,  it.  is  iiiqiortniit  to  nolo  Iwn*  that  the  solutions 
presented  previously  art!  ind«|>oudout  of  the  individual  flow  areas 
at  tin*  inlot  piano.  In  practice,  howovor,  it  is  possible  for  those 
areas,  if  not  pm|M*rly  designed,  to  cause  largo  three-dimensional 
effoots  at  tho  inlot  piano,  due,  for  example,  to  independent  chok- 
ing of  the  primary  stream.  Nevertheless,  it  is  still  either  tho 
nozzle  exit  area  or  geometric  throat,  no/  tkr  inlet  area*,  which 
will  control  the  flow  liehavior.  Thus,  even  if  tliese  three-dimen- 
sional effects  cause  stagnation  pressure  losses  at  the  inlet  plane, 
the  theory  still  gives  the  correct  solution  provided  that  pn  and  pi* 
are  measured  downstream  of  tin*  loss.  In  any  case,  this  situation 
can  lie  avoided  by  designing  the  proper  flow  areas  for  the  inlet 
plane.  This  can  easily  Ik*  done  by  solving  equations  (1)  and  (">) 
simultaneously  at  the  inlet  plane  for  the  flow  parameters  corre- 
sponding to  the  one-dimensional  solution. 

The  procedures  outlined  in  this  section  are  particularly  suited 
to  problems  where  tr-  Vr«M  y/TM  >*  specified.  However,  the 
theory  can  be  applied  equally  as  well  to  generate  solutions  for 
any  suitable  set  of  given  conditions.  The  methods  presented 
here  mav  ulso  1*  extended  to  cases  of  more  than  two  flows. 


Three- Dimniuil  Twe-Streeei  Neule 
Geapeter  Celcelrtiees 


One  of  the  basic  assumptions  of  the  one-dimeusioual  theory  is 
that  the  effects  of  streamline  curvature  are  small  enough  to  lie 
neglected.  This  can  lie  an  important  restriction  in  the  applica- 
tion of  the  tlieory  to  actual  nozzles.  Therefore,  to  demonstrate 
the  effect  of  such  curvatures,  a three-dimensional  two-stream 
nozzle  solution  will  lie  compared  with  the  one-dimensional  theory. 
To  do  this,  an  axially  symmetric  two-flow  nozzle  computation 
was  curried  out.  A brief  description  of  the  computation  is  given 
hen*  to  justify  its  use  in  evaluating  the  one-dimensional  theory. 

The  basic  iden  liehind  this  calculation  can  I rest  lie  understood 
by  recalling  that  the  curves  in  Fig.  2 were  obtained  by  an  integra- 
tion of  equation  (3)  through  tin*  nozzle.  Here  a numerical 
integration  includes  three-dimensional  effects. 

Some  important  assumptions  art*  made  nlmut  the  flow  liefore 
this  solution  is  applied;  namely, 

1 Tlie  system  must  lie  symmetrical  with  respect  to  its  center 
line. 

2 The  two  flows  are  iseutmpic  and  do  not  mix  but  are  in  con- 
tact along  a slipline. 

3 Tlie  static  pressure  is  free  to  vary  along  the  slipline  but 
must  lie  equal  wows  it  at  every  point. 

4 Tlie  primary,  or  inner,  stream  must  lie  everywhere  super- 
sonic while  tlie  secondary,  or  outer,  stream  is  not  restricted. 
13ns  allows  a metliod  of  characteristics  to  lie  n*ed  for  the  primary 
stream  calculations. 


For  any  one  problem,  the  temperature-corrected  mass  flow  ratio, 
the  ratio  of  s|iecifir  heats  for  each  stream,  the  gas  constant  for 
each  stream,  the  coordinates  of  the  wall  shajie,  altd  the  primary* 
stream  inlet  Mach  notnlier  are  roiuodcred  to  Iw  specified. 
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TIm*  flow  field  of  tho  primary  stream  is  constructed  using 
standard  numerical  forms  of  the  method  of  characteristics  equa- 
tions for  axially  symmetric  flow.  The  secondary  stream  is  calcu- 
lated using  one-dimensional,  iscu tropic  flow  relations.  The  flow 
areas  of  the  secondary  stream  am  measured  along  a projection 
normal  to  the  slipline,  Fig.  7.  The  assumption  that  this  flow 
is  one-dimensional  and  has  no  streamline  curvature  effects  is  not 
a serious  restriction.  It  should  lie  noted  that  the  secondary  flow 
passage  occupies  only  a small  (sirtiou  of  the  total  radius  at  any 
section  (relatively  small  mass  flow  rates  and  annular  geometry). 
Therefore,  tlie  streamline  curvature  effects  are  small  compared 
to  those  of  the  primary  stream,  where  they  have  been  accounted 
for. 

The  secondary  and  the  primary  streams  meet  along  the  slipliue. 
The  calculation  of  a succession  of  points  along  this  slipline  involves 
an  iteration  procedure  incorporating  the  calculations  of  the 
flow  properties  of  each  stream  simultaneously,  Fig.  7.  Point  C 
illustrates  a typical  slipline  point  to  lie  calculated.  Point  B,  the 
previous  slipliue  point,  and  point  A,  the  previous  point  on  Mach 
lines  A-R  and  A-C,  are  known.  The  conditions  of  point  C art* 
that  it  be  on  the  Mach  line  A-C  and  that  the  static  pressure  at  C 
lie  the  same  in  the  primary  and  secondary  streams.  These  con- 
ditions are  satisfied  through  a trial-and-error  solution  involving 
the  location  of  point  C. 

As  was  noted,  the  program  is  performing  a numerical  integra- 
tion equivalent  to  the  integration  of  equation  (3).  For  any  set  of 
given  conditions,  the  results  can  lie  presented  in  terms  of  the 
poi/pot  required  to  drive  the  flow.  A reduction  pw/p,,,  corre- 
sponds to  operation  at  a decreased  back  pressure  p. . The  term 
p,  wilt  correspondingly  decrease  until  pw/poi  reaches  a minimum 
value.  Any  value  lielow  this  minimum  will  cause  the  slipline 
iteration  to  converge  on  a secondary  area  ratio,  Aj/.-t,*  < I, 
which  is  physically  meaningless.  This  minimum  value  of  p^/po i 
corresponds  to  choked  flow.  All  values  of  p*»/p* t greater  than 
the  minimum  value  mrre*|xiud  to  uuchoked  flow.  In  short,  the 
three-dimensional  solution  closely  parallels  tlie  nne-dimeusiounl 
solution,  Fig.  4. 

The  three-dimensional  computation  for  Imth  choked  and  mi* 
choked  flows  usually  presents  little  difficulty,  and  tv»  integration 
through  the  geometric  throat  is  generally  smooth.  Difficulties 
may  arise  when  the  accelerating  secondary  stream  approaches 
Mach  1.  This  is  due  to  the  usual  mathematical  sensitivities  en- 
countered in  computing  one-dimeusioual  flow  pm|M*rtir*  near 
Mach  I.  However,  tlie  transition  through  this  region  can  !>* 
made  by  special  technique 

This  study  was  uu«lertakeii  to  find  the  effects  <>f  streamline 
curvature  on  tlie  Is'havior  of  a compound  nozzle.  Tlie  first  cal- 
culation* were  (irrformed  with  geometry  curre.*|M>ndiug  to  tlie 
basic  nozzle  iired  for  tlie  ex|a*rimeutal  portion  of  this  |ia)ier.  Tlie 
flow  fields  in  these  cases  showed  definite  three-dimensional  !*■- 
hnvior,  especially  in  tlie  initial  expansion  region.  However,  tlie 
calculations  »ltowcd  no  important  roaUtring  of  Mark  hnrt,  to  tkat 
tkr  i nrnlropie  <ur<mp/M m wom  ettml tolly  hi/m/. 

To  further  study  three-dimensional  effects,  thr  nozzle  geometry 
was  modified.  A series  of  nozzle  shapes  was  generated  by  apply- 
ing a scale  factor  to  only  the  axial  coordinates  of  the  liasic  nozzle. 
The  resultant  senes  then  consisted  of  nozzles  which  were  identical 
one-dimeimutially  but  had  severe  variation  in  the  slopes  of  the 
contour.  These  geometries  ore  shown  in  Klg.  K. 

The  results  fn»m  this  series  of  computations  are  presented  m 
Ftg.  i*  Curves  a and  h demonstrate  the  compound- choked  flow 
regime.  The  term  p«/pt,  is  plotted  versus  a.**.,  the 

length  from  iht  inlet  l«  the  throat  divided  bv  the  diameter  of  the 
throat.  Notice  that  for  both  curves  the  agreement  with  nm- 
dimenstoiial  results  is  excr'lent  for  relatively  Urge 
ratio*.  For  the  smaller  values  of  the  three-dimm- 

-tonal  effect*  influence  p*  p*,.  This  is  not  surprising  since  in  this 
range  the  nozzle  contour  slope*  are  most  severe  Abo.  the  nozzle 
throat  is  rto«e  enough  to  the  inlet  to  he  in  the  influence  in  the 
initial  expansion  of  the  primary  stream. 
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Curves  c and  d tlemonstratc  the  niiii|*iund-nnrhnkrd  How 
regime.  Again,  agreement  with  the  one-dimensional  results  in 
excellent.  The  discrepancies  f«»r  *lmrt  nozzle*  are  not  so  pro- 
nounced here  since  the  exit  condition*,  am)  not  the  throat  rendi- 
tions, are  dominant.  Some  dejatrture  from  tin*  one-dimensional 
line  in  neen  for  the  longer  nozzles.  It  i*  frit  that  lh>*  may  lie  due 
to  accumulating  inaccuracies  resulting  from  the  finite  difference 
technique*  u*ed.  For  long  nozzles,  tlir»r  error*  may  larrotne 
significant  hy  tlte  time  therxit  plane  i*  reached. 

Tlie  agreement  tielwcen  tlie  onr-dimensionai  and  tltr  axially 
*y  mine  trie  solution*  in  excellent  for  the  r?.ngr  of  nozzle  variation 
examined.  Three-dimensional  effect*  have  little  influence  on  tlie 
level  of  pw/fSt  for  the  nozzle  gcomettie*  considered. 

EiptfiMitai  Rtsitts 

Kxtendve  te*t  jifogram*  conducted  on  a wiile  variety  of  nozzle 
type*  and  geonietne*  have  *hown  excellent  agreement  l«*lween 
the  one-dtmrn*iohal  n»m|*Mind-flow  theory  and  exjiertmenial  re- 
sult*. Thrre  t«wt  program*  have  included  convergent,  rylmdnra), 
and  ronvrrgrtit-divergent  mnilr*  with  U»th  two  ami  three 
stream*. 

The  »urre*M  of  the  one-dimen*Mmal  theory  w not  surprising. 
The  prevaMU  x*rtnm  indurated  that  the  larhavmr  of  rump»Kind 
How  tHiiiin  w muamahly  uiartisiltvr  u>  thrrr~dimm**onai  effect*. 
For  nozzle*  with  »mall  wall  fnrt^uial  effect*,  and  fairly  undis- 
lotted  inlet  flow,  only  the  effect*  of  mixing  ran  rao«e  the  one- 
dunertMUmal  model  to  he  inaccurate  when  applied  to  real  horde*. 
It  t*  reamMialde  to  xentnr  that,  since  the  flow  u titf4*tlent.  muing 
ta  confined  to  a *htar  layer  l «r tween  adjacent  *treatn*  wSirh  grows 
with  axial  pnatina  at  an  angle  ><i  le%*  than  s deg  1^1  - It  *•  there- 
fore rlrar  that  the  shape*  of  thr  .1,  and  the  horde  length  are  U » 
major  factors  which  determine  the  degree  id  mixing.  The  ten- 
dency uf  the  muing  will  lw  to  pump  the  low- velocity  stream*  and 
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to  retard  the  high-velocity  stream*.  The  extent  of  the  influence  of 
mixing  depends  primarily  on  the  flow  rates  of  each  stream.  The 
liehavior  of  streams  with  pro|K>rtionately  low  rates  will  lie  greatly 
affected  hy  mixing  while  those  with  higher  flow  rates  will  cx- 
{lerience  only  small  effects.  Xote,  however,  that  in  thr  important 
roar  of  chokrd  tlou\  mixing  effecta  downatream  of  thr  nozzle  throat 
can  nrrl  no  influence  whatsoever  upon  the  behavior  of  the  flow. 
Tims,  for  many  nozzle  applications,  mixing  will  influence  only  a 
small  |Mirtion  of  the  How. 

licran.se  of  the  space  limitations  of  this  |Hi]ier,  n fully  compre- 
hensive comparison  of  one-dimensional  theory  with  ex|ierimental 
data  front  a wide  variety  of  cont|>ound-How  nozzle  tvjies  and 
geometries  is  inipractieal.  Accordingly,  tlie  basic  nozzle  previ- 
ously discussed  was  selected  for  comparison.  Being  convergent- 
divergent,  it  is  a representative  two-stream  nozzle  in  tliat  it  ran 
display  laith  the  choked  and  uuchoked  regimes  of  roinpoiuid 
flow. 

Test*  were  run  over  a wide  range  »f  U*ni|>eralure-mrrerted 
mass  How  ratios  while  varying  pn/p*  from  approximately  2 to  10. 
'This  allowed  thr  nozzle  to  exhibit  Isith  choked  and  unchoked 
liehavior  at  each  mtw  How  ratio.  Fig.  10  is  a mmjiarisoti  of  ex- 
twrimenlnl  results  with  prediction*  l wised  upon  the  one-dimen- 
sional com]  m mud -flow  tlieory.  llie  rhoketl  flow  regime  i*  the 
straight  |sirtinn  of  tltr  theoretical  line*  ami,  as  rx|iected,  omir* 
at  the  higher  pvt/p *.  Hu*  unrhok«l  H«»w  regime  i*  tlie  curved 
|Mirlinh  of  thr  line*  whirh  orcur*  at  thr  iowrr  primary’  nozzle 
stagnalsm  pressure  ratio*. 

It  ean  la*  *rrn  that  thr  one-dtmrn»tonal  lhr«iry  »how*  rxrrllenl 
agree  men  l with  tlie  rx]irrimri:tal  data,  (airtirularly  in  thr  choked 
flow  regime.  Correlation  with  rxivritnmtal  data  in  the  un- 
rhoktsi  flow  regime  i*  *<>mewhat  lr*»  aemratr.  However,  in  view 
of  tlir  mixing  effect*  previously  di»ru*M*.l,  tht*  was  to  l<*  r\|«retrd. 

I hiring  choked  flow,  all  mixing  downstream  of  the  nozzle  throat 
ran  haw  no  effect  on  thr  flow.  Thu*  the  effective  mixing  length 
for  rhAnl  flow  t»  merely  the  dUtanee  from  the  inlet  |»lane  to  thr 
throat,  tht  the  other  hand,  during  unrhok.'d  flow,  all  mixing 
dowiutream  of  the  twizzle  throat  will  have  a very-  definite  effect 
ott  the  flow  Iwhavmr.  With  unrlnArd  flow,  then,  thr  effective 
muing  length  u the  entire  length  of  the  nozzle. 

A*  an  tin  | tales),  we  alwi  ttlwrrvr  that  ttuxuig  has  little  effect  at 
thr  higher  m.v*»  ffi»w  ratio*  Init  exerts  increasing  influence  a* 
matt*  flow  rata*  ilmvae*.  Furthermore,  the  important  effect 
of  mixing  u to  pump  the  *eoiitdary  flow  ami  therefore  reduce  the 
required  #»*  >« 

It  is  noted  that  the  model  tested  was  not  an  alratlx  designed 
nozzle,  ir,  the  primary  stream  was  independently  rhtAetl  4h| 
•lightly  undrrexpamled  at  thr  inlet  plane.  However,  the  model 
used  in  the  experiment*  was  upnln)  *ufficienliy  i«u  hi  its 
lartilfuptc  dr-  . i rtrttditam*  that  stagnatkin  pressure  l*sse*  in  the 
su  person  it  stream  were  proliahly  not  important,  and  shuck  ksaww 


an*  impossible  in  the  subsonic  secondary  stream.  This  was  indi- 
cated in  the  previous  section,  where  the  three-dimensional  solu- 
tion showed  no  significant  coalescing  of  the  primary  stream  Mach 
lines  nnder  conditions  for  which  the  tests  were  run.  Although 
p.,i /p<n  is  measured  upstream  of  the  inlet,  plane,  these  argu- 
ments justify  our  neglect  of  total  pressure  losses  in  all  calcula- 
tions. 

Concluding  Remarks 

A new  one-dimensional  llienry  describing  the  liehavior  of  cont- 
jioii  nil-compressible  nozzle  flows  hits  Iwen  developed  and  its 
implications  have  lieen  examined  from  a munlicr  of  viewpoints. 

The  theory  yields  simple  algebraic  methods  for  calculating  the 
o|«> ration  of  eonqiouiid-compressible  nozzles.  Comparison  of  the 
algeluaic  results  with  those  of  three-dimensional  flow  field  com- 
putations indicates  that  the  effects  of  streamline  curvntnre  are 
not  iiu|Mirtant  for  many  practical  nozzle  configurations.  Com- 
parison of  the  ulgehraic  results  with  ex|>erinientul  data  for  flows 
with  i»iim|M>rlant  mixing  effects  shows  that  the  theory  can 


accurately  predict  the  liehavior  of  real  devices. 

A complete  definition  of  tlie  range  of  applicability  of  the  simple 
theory  requires  a great  deal  of  ex|ierimeutal  experience.  For  ex- 
ample, little  information  is  available  alniut  the  rate  of  growth  of 
the  mixing  zone  lietwcen  the  high-velocity  streams,  and  no  in- 
formation is  available  about  the  detailed  nature  of  compound- 
shock  waves.  Furthermore,  no  predictions  about  the  thrust  or 
nozzle  efliriency  can  be  made  when  the  flow  is  choked  until  the 
losses  caused  by  a compound-shock  wave  are  known.  Conse- 
quently, it  appears  that  the  compound-compressible  nozzle  pro- 
vides a iiumtier  of  interesting  and  important  areas  of  research. 
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